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^ • Abstract: We confront a very wide body of HERA diffractive electroproduction data with 

^ • the predictions of the colour dipole model. We focus upon three different parameterisations 

D , of the dipole scattering cross-section and find good agreement for all observables. There can 



be no doubting the success of the dipole scattering approach and more precise observations 
are needed in order to expose its limitations. 
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1. Introduction 

In the colour dipole model |l], ^, the forward amplitude for virtual Compton scattering 
is assumed to be dominated by the mechanism illustrated in Fig.|| in which the photon 
fluctuates into a qq pair of fixed transverse separation r and the quark carries a fraction z 
of the incoming photon light-cone energy. Using the Optical Theorem, this leads to 

a!;il = J dzd\ \^^''^[r,z,Q^)\''a{s\r) (1.1) 

for the total virtual photon-proton cross-section, where are the appropriate spin- 

averaged light-cone wavefunctions of the photon and a{s*,r) is the dipole cross-section. 
The dipole cross-section is usually assumed to be independent of z, and is parameterised 
in terms of an energy variable s* which depends upon the model. 

Thus using Eq.(|l.lD we are able to compute the deep inelastic structure function 
F2{x,Q'^). The power of the dipole model formulation lies in the fact that the same dipole 
cross-section a{s* ,r) appears in a variety of other observables which involve the scattering 
of a real or virtual photon off a hadronic (or nuclear) target at high centre-of-mass (CM) 
energy. The largeness of the CM energy guarantees the factorization of scattering ampli- 
tudes into a product of wavefunctions and a universal dipole cross-section. In this paper we 
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Figure 1: The colour dipole model for elastic Compton scattering — > ^*p. 

wish to test the universahty of the dipole cross-section using a wide range of high quality 
data collected at the HERA ep collider. Moreover, we also wish to examine the extent 
to which to data are able to inform us of the role, if any, played by non-linear saturation 
dynamics. 

Specifically, we shall consider three particular parameterisations of the dipole cross- 
section which have been presented in the literature and which will be discussed in more 
detail below. All three have been constrained by fitting only to the HERA data on F2{x, Q^) 
and so they can be used to make genuine predictions for other observables. In this paper 
we will compare those predictions to data on the charm structure function F2{x,Q^), 
the cross-section for Deeply Virtual Compton Scattering (DVCS), the cross-section for 
diffractive J/"^ production and the diffractive structure function F^^^\ 

In the first instance our results confirm the validity of the dipole model: it does make 
sense to speak of a universal dipole cross-section which is able to account for both soft 
and hard diffraction in a wide variety of photo-processes. We also attempt to ascertain the 
extent to which the new data are able to discriminate between the predictions of the three 
dipole models we use. These models can perhaps best be described as parameterisations 
which incorporate certain general theoretical ideas and can be viewed as providing an indi- 
cation of the uncertainties remaining in the dipole cross-section once the precise structure 
function data have been accounted for. Unfortunately we shall find that the new data are 
not quite precise enough to discriminate between the models or to add any significant evi- 
dence for saturation beyond that already present in the F2 data. We do however find that 
in diffractive photo/electro-production progress is hindered by the lack of a precise enough 
measurement of the forward slope parameter B which determines the t-dcpendence of the 
final state proton. We take this parameter from data and its error translates into an uncer- 
tainty on the normalisation of the predicted cross-sections. A more precise measurement 
of this quantity would provide a significant additional constraint. 

It is clear that having extracted the dipole cross-section from data, one would hope 
eventually to explain the result using QCD. Unfortunately, while QCD in its present state 
of development is able to suggest qualitative features of the dipole cross-section, more 
quantitative predictions are not possible without severe approximations. To remedy this 
situation is a challenge which lies beyond the scope of this paper. 
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2. The dipole cross-section 



We now turn to the three different models used to describe the dipole cross-section. Before 
doing so however, we shall first discuss our choice of photon wavefunction. For small r, the 
light-cone photon wavefunctions are given by the tree level QED expressions Q 

\^^{r,z,Q^)\^ = -^a,^^e2Q2^2(^ _ ,fKl{er) (2.1) 
^ /=i 



3 ^ 



^^ir,z,Q')\' = ^«en.E4 {[^' + (1 - zf]e'K!{er)+m}K',{er)} (2.2) 

f=i 



27r2 

where 



= z{l - z)Q^ + mj . (2.3) 

Here Kq(x) and Ki{x) = —dxKQ{x) are modified Bessel functions and the sum is over 
all = 4 quark flavours /. These wavefunctions decay exponentially at large r, with 
typical r-values of order at large and of order "rnj^ at = 0. However for 
large dipoles r > 1 fm, which are important at low Q^, a perturbative treatment is not 
really appropriate. In this region some authors modify the perturbative wavefunction 
by an enhancement factor motivated by generalised vector dominance (GVD) ideas ||^, 
m, while others 0] achieve a similar but broader enhancement by varying the quark 
mass^. In practice the difference between these two approaches only becomes important 



when analysing the precise real photoabsorption data from fixed-target experiments [10|. 
Since we will not consider these data here, we will adopt the simpler practise of using a 
perturbative wavefunction at all r-values, and adjusting the quark mass to fit the data. 

Turning now to the dipole cross-section, all three models are consistent with the 
physics of colour transparency for small dipoles and exhibit soft hadronic behaviour for 
large dipoles. As stated above, the model parameters are determined by fitting only to 
the DIS structure function data. The resulting dipole cross-sections can then be used to 
make genuine predictions for other reactions. Since the details of all three models have 
been published elsewhere, we shall here summarise their properties only rather briefly. 

2.1 The FS04 Regge model 

This simple model Q combines colour transparency for small dipoles r < tq with "soft 
pomeron" behaviour for large dipoles r > ri by assuming 

a{xm,r) = Aur'^x^" for r < ro and 

= Asx:;:^^ for r > n, (2.4) 



where 



q2 f ^ ^ 4rnj 



Q2^W2 \ Q 



1 + ^ . (2.5) 



^For a fuller discussion of these points see 
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For light quark dipoles, the quark mass my is a parameter in the fit, whilst for charm quark 
dipoles the mass is fixed at 1.4 GeV. In the intermediate region tq < r < ri, the dipole 
cross-section is given by interpolating linearly between the two forms of Eq.(p.4|). 

If the boundary parameters tq and ri are kept constant then this parameterisation 
reduces to a sum of two powers, as might be predicted in a two pomeron approach, and 
can be thought of as an update of the original FKS Regge model Q to accommodate the 
latest data. It is plainly unsaturated, in that the dipole cross-section obtained at small 
r-values and fixed grows rapidly with increasing s (or equivalently with decreasing x) 
without damping of any kind. 

2.2 The FS04 Saturation model 

Saturation can be introduced into the above model by adopting a method previously utilised 
in |T^]. Instead of taking tq to be constant, it is fixed to be the value at which the hard 
component is some specified fraction of the soft component, i.e. 

cr{xm,ro)/(7{xm,ri) = f (2.6) 

and / instead of is treated as a parameter in the fit. This introduces no new parameters 
compared to the Regge model. However, the scale rg now moves to lower values as x 
decreases, and the rapid growth of the dipole cross-section at a fixed, small value of r begins 
to be damped as soon as tq becomes smaller than r. In this sense we model saturation, 
albeit crudely, with ro the saturation radius. 

2.3 The CGC saturation model 

In addition we shall consider the CGC dipole model originally presented by lancu, Itakura 
and Munier [|l^]. This model aims to include the main features of the "Colour Glass 
Condensate" regime, and can be thought of as a more sophisticated version of the original 
"Saturation Model" of Golec-Biernat and Wiisthoff Q. Since the original lancu et al dipole 
cross-section was obtained after a three fiavour fit to the DIS data it is not well suited to 
making predictions for processes involving charm quarks. Consequently, we instead use a 
new four-fiavour CGC fit due to Kowalski, Motyka and Watt |l^] . 



The parameters of the FS04 models were determined by fitting the recent ZEUS F2 
data 1 14] in the kinematic range 



0.045 GeV^ < < 45 GeV^ x < 0.01 



(2.7) 



whilst the CGC fit of pi is to data with > 0.25 GeV^ (the other hmits are as for 



FS04). The corresponding HI data |15| could also be used, but it would then be necessary 
to fioat the relative normalisation of the two data sets. We do not do this since the ZEUS 
data alone suffice. The resulting parameter values are tabulated in the original papers; we 
do not reproduce them here, but confine ourselves to some general comments. 
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Both the FS04 saturation and the CGC models gave excellent fits to the F2 data, 
while the FS04 Regge fit was not satisfactory, suggesting that saturation may be required 
to fit the data However this evidence for saturation depends upon using low-Q^ data 
and disappears if the data are restricted to > 2 GeV^, whereupon excellent fits can be 
obtained in all three models. 

We use the CGC fit with do = 35.7 mb presented in Table 5 of |13| and note that 



although the fit is to data with > 0.25 GeV^ the fit is actually very good all the way 
down to Q2 = 0.045 GeV^. 

At this point we have three well-determined parameterisations of the colour dipole 
cross-section. These can be used to yield predictions for other processes. In the next 
sections we shall take a look at Deeply Virtual Compton Scattering (DVCS), the charm 
structure function exclusive J/\I' production and the diffractive structure function 
F^^^\ We always choose to show the Regge fit, even though it does not fit the F2 data 
particularly well, in order to indicate the discriminatory power of the data. We stress 
that in all cases, the photon wavefunctions and dipole cross-sections are precisely those 
determined from the fits to F2 data, without any adjustment of parameters. 

3. Exclusive processes 




Figure 2: The colour dipole model for the exclusive reactions 7*p Ap. 

We are interested in the exclusive processes 

Y+P^A + p A = 'y,p,J/'^,.... (3.1) 

In the dipole model they occur via the mechanism of Figj2| and are described by amplitudes 
which satisfy 

ImAxis, t = 0)=s^j dzd\ ^^-^{r, z)* ^J^(r, z, Q^) <j{s\r) (3.2) 

hji 

where X = L,T for longitudinal and transverse photons respectively, h and h are the quark 
and antiquark helicities and ^'^^(r, z) is the light-cone wavefunction of the particle A. The 

^The FS04 saturation fit has a = 155 for the 156 data points we consider. For the same data, the 
CGC model has = 160 and the FS04 Regge fit has x^ = 428. 
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forward differential cross-section is then given by 



da 
H 



^ \A^{s,t = 0)\\l+P^) , (3.3) 



where the correction from the real part of the amplitude 

ReAx{s,t = 0) 



linAxis,t = 0) 



can be estimated using dispersion techniques. Predictions for the measured total cross- 
sections are then obtained by assuming an exponential t-distribution and integrating over 
t to obtain 



B dt 



(3.4) 

i=0 



where the value of the slope parameter B is taken from experiment. We refer to |13, 16 
for models which attempt a more sophisticated treatment of the dependence upon the 
momentum transfer, t. 

The exclusive processes we consider in this paper are DVCS and J/^ production. Light 



meson production was studied in our previous paper ||17|] where we found that the dipole 
model predictions generally agree well with the data modulo a rather strong dependence 
upon the meson wavefunction. For an alternative investigation of the link which exists 
between low x DIS and exclusive processes at high energies we refer to 

3.1 Deeply virtual Compton scattering 

In deeply virtual Compton scattering 

7*+P^7 + P, (3.5) 

the final state particle is a real photon and dipole models provide predictions for the 
imaginary part of the forward amplitude with no adjustable parameters beyond those used 
to describe DIS. To calculate the forward cross-section a correction for the contribution of 



the real part of the amplitude has to be included. This correction was estimated in [19| 
and found to be less than « 10% and of a similar size in different dipole models. Here 
we shall estimate the correction using the FS04 Regge model, where the real part is given 
by the Regge signature factors. Predictions for the measured total cross-sections are then 
obtained using Eq.(|3.4]) where the value of the slope parameter B is taken from experiment. 



The predictions of all three models are compared with the HI data |2C] in Fig.H and 



Fig.| and with the ZEUS data ^ in Fig.| and Fig.|.3 For the HI data we use a fixed 
value B = 6.02 GeV^^ for the slope parameter which is in accord with the HI measurement 
oi B = 6.02 ± 0.35 ± 0.39. For the ZEUS data we take B = A GeV~2 which is compatible 
with their data. Bearing in mind this normalisation uncertainty, the agreement is good 
for all three models, although significant differences between the models appear when the 
predictions are extrapolated to high enough energies, as one would expect. 

^Note that throughout this paper the curves labelled 'FS04 no sat' correspond to the predictions of the 
FS04 Regge model. 
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Figure 3: Comparison of tiie HI DVCS data ||2^ witli tlie predictions of the three models discussed 
in the text: dependence at W = 82 GeV. 
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Figure 4: Comparison of the HI DVCS data ||2^ with the predictions of the three models discussed 
in the text: W dependence at = 8.0 GeV^. 

4. J/^/ and the charm structure function 

We now move on to the predictions for the charm structure function and for exclusive J j'^ 
production 

7* +p J/^- +p . (4.1) 

All three models assume a charm mass of mc = 1.4 GeV when fitting the F2 data. Since 
the exclusive process is rather more sensitive to the charm mass we will allow mc to vary 
a little without adjusting the dipole cross-section. This is permissible for small enough 
variations. 

4.1 The charm structure function 

We begin by discussing the charm structure function ^2(^,(5^), since the results are in- 
dependent of the vexed question of the vector meson wavefunction. The charm structure 
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Figure 5: Comparison of tiie ZEUS DVCS data with the predictions of the three models 
discussed in the text: dependence at W = 89 GeV. 
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Figure 6: Comparison of the ZEUS DVCS data with the predictions of the three models 
discussed in the text: W dependence at = 9.6 GeV^. 



function is given by 



4^2 



ar. 



/ L T \ 



where in calculating the total virtual photon-proton cross-sections using Eg. ( |1 . ![ ) , only 
the charm component of the light-cone wavefunctions ( |2.1| ) is retained. The resulting 
predictions are compared to the ZEUS [^] and HI [23| data in Figs.^-^. A good account 
of the charm structure function data can be obtained in all three models by choosing values 
for the charmed quark mass in the reasonable range 1.3 < rUc < 1.5 GeV. The key question 
is whether one can obtain similarly accurate predictions for the J/^ electroproduction data 
using charm mass values in the same range. 
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Figure 7: Comparison of the FS04 saturation model predictions for the charmed structure function 
with data ^ ||. 



4.2 J/^* wavefunctions 

To calculate vector-meson production we need to know the light-cone wavefunctions of the 
vector mesons. Three different types of vector-meson wavefunction were studied for light 



mesons in [17|. For the J/^ we shall use wavefunctions of exactly the same functional 
form, but with the parameters adjusted to take into account the mass and charge of the 
charmed quark and the experimental value of electronic decay width of the J/^'-meson. In 
what follows, we shall comment briefly on the resulting wavefunctions, referring to jl^] for 
detailed formulae and discussion. 



In the DGKP approach |24], the r and z dependence of the scalar wavefunction is 



assumed to factorise into a product of gaussians. In the other two cases considered in [17|, 
it is obtained by taking a given wavefunction in the meson rest frame. This is then boosted 
into a light-cone wavefunction using the Brodsky-Huang-Lepage prescription, in which the 
expressions for the off-shellness in the centre-of-mass and light-cone frames are equated 
p5| ] (or equivalently, the expressions for the invariant mass of the qq pair in the centre- 
of-mass and light-cone frames are equated ^^). In the simplest version of this approach, 
the wavefunction assumes a gaussian form in the meson rest frame. Alternatively NNPZ 



1 27] have supplemented this by adding a hard "Coulomb" contribution in the hope of 



improving the description of the rest-frame wavefunction at small r. However there are 
theoretical problems with this latter wavefunction, as discussed in |17], and we shall confine 
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Q"=12.0 Q"=18.0 Q" = 25.0 = 30.0 




our discussion here to the DGKP and "Boosted Gaussian" wavefunctions. 

The appropriate wavefunction parameters for the DGKP and Boosted Gaussian wave- 
functions in the J/^' case are given, in the notation of [^], in Tables |l] and ^ respectively.^ 
For a charm quark of given mass, they are chosen so that the wavefunction is normalised 
and so that the predictions of the decay constant are consistent with the experimental 
value fj/^ 



0.273 ± 0.005 GeV (see Q for further details). 

DGKP parameters 







ujt 


Ml 


Nt 


1.4 


0.688 


0.560 


18.941 


8.280 


1.35 


0.688 


0.568 


18.941 


8.616 



Table 1: Parameters and normalisations of the DGKP light-cone wavefunctions in appropriate 
GeV based units. 

The resulting wavefunctions are shown in Figs.|lO[-[lT| for the case rric = 1-4 GeV. Like 
the corresponding wavefunctions for the p and (/> mesons (shown in [|l^ ) , the wavefunctions 
peak at z = 0.5 and r = 0, and go to zero as z — > 0, 1 and r ^ oo. As expected, for the 
J/^ case the peaks are much sharper. We also see that the DGKP and Boosted Gaussian 

^Compare with Tables 2, 3 and 4 of 0. 
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Figure 9: Comparison of the CGC model predictions for the charmed structure function F2 with 





Figure 10: The J/^I'-wavefunctions I^J^^p (left) and (right) in the DGKP model. 



wavefunctions are qualitatively similar, with the transverse wavefunction having a broader 
distribution than the longitudinal wavefunction in each case. 



11 



Figure 11: The J/^'-wavefunctions j^f^p (left) and j'l'-^p (right) in the Boosted Gaussian model. 

Boosted Gaussian parameters 







Ml 


Mt 


ML) 


fv{T) 


1.4 


2.44 


0.0363 


0.0365 


0.262 


0.293 


1.35 


2.44 


0.0369 


0.0370 


0.266 


0.288 



Table 2: Parameters and normalisations of the Boosted Gaussian light-cone wavefunctions in 
appropriate GeV based units. 



4.3 J/^* electroproduction 

Given the forms of tiie photon and J/^' wavefunctions, the cross-section in Eq.(|3.4[) can be 



calculated from Eqs.(3.2)~(p.4|). However, to do this requires an estimate of the correction 
due to the real part in Eq.( |3.3|) and of the slope parameter B in Eq.( |3.4[ ). 

The correction factor (1 + /3^) for the real part of the amplitude is, as for the DVCS 



case, determined from the FS04 Regge model. The resulting values are illustrated in Fig. 12 
which shows the values obtained as functions of W at = 0, 10 GeV^. As can be seen, 
the corrections from the real parts in Eq. (|3.3| ) are a significant, but not large, correction. 



For the slope parameter B in Eq.(3^) we use the simple parameterisation (in GeV units) 



B = nI lilL^ + i) (4.2) 

with N = 0.55 GeV"2 for the J/^' which is in accord with the data [M, ^9(1 . 

Having fixed the real parts and the slope parameter, predictions for the production 
cross-section can be made for both the Boosted Gaussian and DGKP wavefunctions. The 
charmed quark mass is adjusted within the range allowed by the charm structure function 
data, i.e. we use rric = 1.4 GeV for the FS04 Regge model and nic = 1-35 GeV for the two 
saturation models. 
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Figure 12: The W dependence of the squared ratio of the real to imaginary parts of the J/'i' 
production amplitude (off transverse photons) obtained using the FS04 Regge dipole model. 



The predictions for the dependence of the total cross-section (jTOT _ _|_ ^^l 
take e = 0.98) are presented in Figjl3|. There is very good agreement between theory and 
data with relatively little dependence upon the choice of meson wavefunction.^ We note 
that for ^5 GeV^, the predictions are very sensitive to the choice of charm quark mass. 
For example, the small difference between the predictions of the CGC and the other two 
model predictions in this region can be eliminated by fine-tuning the charm quark mass in 
the CGC case from 1.35 to 1.32 GeV. 

The ly-dependence is shown for a range of values in Figs. 14 -p^. As before, agree- 
ment is good and there is little dependence on the choice of meson wavefunction. Large 
differences between the models only arise at energies beyond the current experimental range 
{W > 400 GeV), but again the differences between the two saturation models can be sig- 
nificantly reduced by fine-tuning the chosen values for the charmed quark mass, especially 
at low Q^. 

Finally, our predictions for the dependence of the cross-section ratio R = cjl/'^t 
at = 90 GeV are shown in Fig. 16. There is once again good agreement between theory 



and data with slightly more dependence upon the choice of meson wavefunction. Note 
that the longitudinal-to-transverse ratio R increases approximately linearly with and, 
in contrast to the p and (j) cases [^], does not flatten out in the currently accessible range 
of Q^. This is to be expected, since in the extreme non-relativistic limit the wavefunction 
approaches a delta function whence the ratio R oc Q"^ /ml for all . 



^Unlike the case for the hght mesons 
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Figure 13: Comparison of the model predictions to the data ^ for exclusive J/"^ meson 
production: dependence. Upper plot: DGKP wavefunction. Lower plot: Boosted Gaussian 
wavefunction. 
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Figure 14: Comparison of tlie model pre(iictions to the <iata |g8|, ^ for exclusive J meson 
proiiuction using the DGKP meson wavefunction: W depen(ience. 
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Figure 15: Comparison of tlie model predictions to the data ^ for exclusive J/^f meson 
production using the Boosted Gaussian meson wavefunction: W dependence. 
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Figure 16: Comparison of tlie model predictions to the data |2^, ^ on the ratio R for exclusive 
J/vf meson. 
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5. Diffractive deep inelastic scattering (DDIS) 



l-z -k, 
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Figure 17: The qq and qqg contributions to F^''^\ 

To conclude our study we turn to the diffractive deep inelastic scattering (DDIS) 
process 

J* + p ^ X + p , 

where the hadronic state X is separated from the proton by a rapidity gap. In this process, 
in addition to the usual variables x and there is a third variable M^. In practice, x 
and M|- are often replaced by the variables xp and /?: 

+ ^-^"^M^ +Q2 • ^^-^^ 

In the diffractive limit s ^ Q^, m\ and so xjp ^ 1 . 

In the dipole model, the contribution due to quark-antiquark dipoles to the structure 
function F^^^^ can be obtained from a momentum space treatment as described in pO| , 0. 
However, if we are to confront the data at low values of (3, corresponding to large invariant 
masses Mx, it is necessary also to include a contribution from the higher Fock state qqg. 
We can estimate this contribution using an effective "two-gluon dipole" approximation due 



to Wiisthoff pg], as illustrated in Fig, 17 



Again, the predictions obtained in this way involve no adjustment of the dipole cross- 
sections and photon wavefunctions used to describe the F2 data. We are however free 
to adjust the forward slope for inclusive diffraction (B) within the range acceptable to 
experiment, which means that the overall normalisation, but not the energy dependence, 
of F^^^^ is free to vary somewhat. We take B = 6.8 GeV~^ when making our CGC and 
FS04 saturation predictions and B = 8.0 GeV~^ when making the FS04 Regge predictions. 
Note that a value of 8.0 GeV~^ is rather high compared to the ~ 6 GeV~^ favoured by 
the HI FPS data |l| although it is in the range allowed by the ZEUS LPS data |3|. The 
need for a larger value of B for the FS04 Regge model arises since the corresponding dipole 
cross-section is significantly larger than the FS04 saturation model at large values of r and 
this enhancment is magnified in inclusive diffraction since it is sensitive to the square of 
the dipole cross-section. We should also bear in mind that the tagged proton data are 
subject to an overall ~ 10% normalisation uncertainty. We are also somewhat free to vary 
the value of as used to define the normalisation of the model dependent qqg component, 
which is important at low values of p. Rather arbitrarily we take Os = 0.1 and take the 
view that the theory curves are less certain in the low /? region. 
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Figure 18: Model predictions compared to the ZEUS LPS data 



In FigJT^ we compare the recent ZEUS LPS data [^] on the xp dependence of the 
structure function F^^^^ at various fixed and /? with the models^. The agreement is 
good except at the larger xjp values. Indeed, the values per data point are very close 
to unity for all three models for xjp < 0.01. Disagreement at larger xjp is to be expected 
since this is the region where we anticipate a significant non-diffractive contribution which 
is absent in the dipole model prediction. Note that the three models produce similar 
predictions at larger values of /?. 



Contamination from secondary exchanges is avoided in the FPC data |34], in which 
the non-diffractive contribution is explicitly removed by the ZEUS method of analysis. 
However this comes at the cost of including proton dissociation contributions, since the 
mass of the target fragments (My) is only limited to being below 2.3 GeV. The predictions 



of our three models are compared with these data in Fig.|l^ and Fig. 20, where the theory 
curves have been divided by 0.7 to allow for proton dissociation contributions, assumed to 
be a fixed fraction of the total, independent of the kinematic variables. This is obviously a 
crude approximation but it is supported by the experimental analyses. These data suggest 
that the CGC and Regge model predictions are perhaps overshooting the data at low /3 
whilst the three models produce very similar predictions at larger values of /?, i.e. /? > 0.4. 



'Predictions for the original lancu et al CGC model have previously been published in 
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Figure 19: Model predictions compared to the ZEUS FPC data with My < 2.3 GeV (low Q^) 



Specifically, for the FS04 saturation model = 238 for the entire 156 data points "whilst 
the CGC and FS04 Regge models have = 147 and 150 respectively for the 94 data points 
with (3 > 0.3. Ho"wever, "we must be careful not to overinterpret the data. The low P region 
is precisely the region where there is an appreciable qqg component which is subject to a 
rather large theoretical uncertainty. Moreover, we must not forget that (for the FPC data) 
we have assumed that the fraction of events which contain a dissociated proton is constant. 
Recall also that there is an uncertainty in the value of the forward slope parameter used 
in determining the theoretical predictions. 



Comparison to the HI data with tagged protons [31 1 is to be found in Fig. 21 and 
Fig.^ The story is similar to that for the ZEUS data and the evidence for an overshoot of 
the CGC and Regge model predictions at low (3 is strengthened. The = 2.7 GeV^ panes 



in Fig, 22 illustrate this point the best. Again, we should not interpret this as evidence 
against these dipole models due to the uncertainty in the qqg contribution in the low /3 
region. The agreement between all models and the data at larger values of f3 and low 



-20- 



■ ZEUSFPC 

FS04 sat b = 6.8 GeV 




10"^ 10"^ 10"^ 10"^ 10"' 10"^ 10'^ 10"' 10"^ 10"" 10"' 10"^ 10"^ 



Q' = 8GeV' Q' = 14GeV' Q' = 27 GeV' Q' = 55 GeV' 



Figure 20: Model predictions compared to tiic ZEUS FPC data witli My < 2.3 GeV (high Q^) 



enough xjp is satisfactory.^ 

In Figsj2^-p^ we compare to the HI data coUected with some proton dissociation, 
i.e. My < 1.6 GeV [^], and in order to include the effect of proton dissociation we now 
divide the theory by a factor 0.8. The FS04 saturation model performs well over the whole 
kinematic region with xp < 10^^ whilst once again the CGC and Regge models overshoot 
the data at low /?. On these plots we also show the result of computing F^^^^ without 
including the qqg component: clearly the data require an additional contribution at low 
p. We note that it is possible to improve the quality of the agreement between data and 
theory if we are allowed to increase the fraction of proton dissociation assumed in the data 
(there is perhaps a hint that this increase should also be slightly larger at higher values 
of Q^). For example, for the 170 data points in Fig,2(: with /3 > 0.4 and xp < 0.01 we 
find values of 244, 330 and 263 for the FS04 saturation, FS04 Regge and CGC models 
respectively if we globally decrease the scaling factor from 0.8 to 0.7. 



■^The FS04 saturation fit has = 37 for the 40 points with xjp < 0.01 whilst the CGC and FS04 Regge 
fits both have P^r data point below 2 for the 18 points with /3 > 0.35. 
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Figure 21: Model predictions compared to the HI FPS data: xp dependence pn. 



Before leaving the HI data we should say that, strictly speaking, all of the HI data 



refer to the reduced diffractive cross-section, a. 



D(3) 



quantity and the diffractive structure function F2 
than 10%. 



D(3) 



However the difference between that 
is mostly negligible and is never more 



In summary, the DDIS data at large enough /5 > 0.4 and small enough xp < 0.01 are 
consistent with the predictions of all three dipole models. However the data themselves 
would have a much greater power to discriminate between models if the forward slope 
parameter were measured to better accuracy. At smaller values of /?, the data clearly reveal 
the presence of higher mass diffractive states which can be estimated via the inclusion of 
a qqg component in the dipole model calculation under the assumption that the three- 
parton system interacts as a single dipole according to the universal dipole cross-section. 
The theoretical calculation at low (3 must be improved before the data in the region can be 
utilised to disentangle the physics of the dipole cross-section. Nevertheless, it is re-assuring 
to observe the broad agreement between theory and data in the low f3 region. 
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Figure 22: Model predictions compared to the HI FPS data: f3 dependence 
6. Conclusion 

The dipole scattering approach, when apphed to difFractive electroproduction processes, 
clearly works very well indeed. The HERA data now constitute a large body of data which 
is typically accurate to the 10% level or better, and without exception the dipole model 
is able to explain the data in terms of a single universal dipole scattering cross-section. 
Perhaps the most important question to ask of the data is the extent to which saturation 
dynamics is present. Although the F2 data suggest the presence of saturation dynamics 
the remaining data on exclusive processes and on F^^^^ are unable to distinguish between 
the models we consider here: these data are therefore unable to offer additional informa- 
tion on the possible role of saturation. We do note that a more accurate determination of 
the forward slope parameter in diffr active photo/electro-production processes would sig- 
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Figure 23: Model predictions compared to the HI data with My < 1-6 GeV (3 dependence 
at Xip = 0.0003. Green dotted curve shows the contribution without including the qqg component. 

nificantly enhance the impact of the data. However, it is hard to avoid the conclusion that 
only with more precise data or with data out to larger values of the centre-of-mass energy 
will we have the chance to make a definitive statement on the role of saturation without 
the inclusion of the low F2{x, Q^) data in the analysis. 
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Figure 26: Model predictions compared to the HI data with My < 1-6 GeV [35]: xjp dependence. 



